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spherical capsules with the latent heat storage system in east-west oriented greenhouse, is analyzed and 
discussed. The excess heat in the greenhouse was stored in the packed bed through the diurnal period and 
extracted at night. An experimental comparative study was conducted in two greenhouses installed in the 
Research and Technologies Centre of Energy (CRTEn) in Tunisia. The greenhouse heat balance in different 
Keywords: components of the greenhouse with a heating system (cover, canopy, soil, inside air and packed bed solar 
Latent heat storage collector), has been used to investigate the impact of the phase change material (PCM) on the greenhouse 
Packed bed sölar-air heaters temperature and humidity. Results show that the nocturnal temperature inside the greenhouse equipped 
Control climate greenhouse : a : A : 

by a heating system exceeds the temperature inside the conventional greenhouse by 5 °C. The relative 
humidity was found to be of an average 10-20% lower at night time inside the heated greenhouse. 
The nighttime recovered heat of the solar system attained 31% of the total requirements of heating. An 
economic evaluation was made to find the payback period of the new solar air heater, the potential in 
saving fuel energy and the reducing carbon dioxide emissions. 
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1. Introduction 


Greenhouses have been used in various forms of protection for 
the plants when it would have been impossible to grow outside 
because of the adverse climate, diseases and pests. The greenhouse 
permits to achieve maximum returns from cultivation by creating 
better growing conditions. The growth factors namely light, humid- 
ity, temperature and nutrition should be delivered and maintained at 
optimal levels [1,2]. To keep the optimum growing environment, 
greenhouses should allow high light transmittance, sufficient venti- 
lation efficiency, low construction costs and low heat consumption 
[1,3]. With the present increase in the cost of fossil fuels, alternate 
energy sources may help to increase the profitability of greenhouses 
and the income of farmers [4]. The basic strategy of the greenhouse 
passive heating system is to reduce the heat losses while transferring 
the excess heat from inside the greenhouse during the day to heat 
storage. This heat is used during the night to satisfy the heating 
needs of the greenhouse [5,6]. 

Alternate sources of energy like solar energy have been 
proposed for heating greenhouses and providing the optimum 
inside conditions during winter months. Solar energy is a uni- 
versally available source, but its practical use still presents some 
technical, and most importantly, economic problems [7]. Some 
solar systems have been developed but they are too expensive for 
commercial use [7]. Several types of passive solar systems and 
techniques have been proposed and used by [8-11]. Bargach et al. 
[10] used a solar fat-plate collector to improve the internal 
greenhouse microclimate and tested it under two types of days; 
a clear day and a cloudy one. Du et al. [12] developed a model to 
predict the performance of the greenhouse that is heated with a 
heat-pipe system. The model is validated with experimental data. 

Thermal energy can be stored as sensible heat, latent heat, reaction 
heat or combination of those forms. Many thermal energy storage 
systems and applications are cited in literature [13-17]. The most 
important existing greenhouse heating systems that use a sensible 
stored energy are: water storage [18], rock bed storage [19,20] and 
ground storage [21,22]. The latent thermal storage has many advan- 
tages over the sensible one such as: high heat capacity, less volume, 
low storage temperature, thermal energy stored and released at an 
almost constant temperature. There is a large number of Phase Change 
Materials (PCMs) that melt and solidify at a wide range of tempera- 
tures which are used in greenhouse, such as CaCl- 6H20, 
NazSO,- 10H20 and paraffin [23-27]. Kern and Aldrich [28] utilized 
a CaCl -6H20 in aerosol cans that was used to investigate energy 
storage possibilities both inside and outside a 36 m? ground area of 
the greenhouse. The energy storage unit inside the greenhouse 
(contained 1650 kg PCM) absorbed the excess energy of warm air 
inside the greenhouse during the daytime. The energy storage unit 
outside the greenhouse contained 1376.4 kg PCM and two solar air 
collectors with 8.55 m? surface area each. Results of this study showed 
that, the energy stored by the outside unit was between 105.5 and 
158.25 MJ, these were, 21.1 and 31.65 MJ for the internal unit. Huang 
and Toksoy [29] placed 3000 kg of CaClz- 6H20 in the heat exchanger 
inside a 200 m? glass covered greenhouse, and found that the interior 
air temperature was maintained at 2 °C higher than the outdoor 
conditions. Levav and Zamir [30] used CaCl - 6H20 in greenhouses and 


reported that the required air temperature in the greenhouse was 
achieved without any increase in the relative humidity. Boulard and 
Baille [31] utilized, in a greenhouse with 176 m? ground areas, double 
polycarbonate-cover and forced ventilation, 2970 kg of CaCl,6 H20. 
The PCM was packed in containers and placed along the north wall. 
The PCM could provide 30% night heating needs during the winter 
period. Boulard et al. [32] used a CaCl,-6H20 latent heat storage 
system in a double-skin polycarbonate greenhouse equipped with a 
forced ventilation system. They found the inside greenhouse tempera- 
ture roughly 10 °C higher than the outside temperature during typical 
nights of March and April. Ozturk et al. [33] performed an experi- 
mental evaluation of energy and exergy efficiency of a seasonal latent 
heat storage system for greenhouse heating using 6000 kg of paraffin 
wax as a PCM in a 180 m? greenhouse ground areas. They observed 
that the average energy and exergy efficiency of the arrangement was 
40.4% and 4.2%, respectively. 

The aim of this work is to determine the provided heat by a 
Solar Air Heater with Latent Storage Collector (SAHLSC) inside an 
Insulated Greenhouse with Latent Heat System (IGLHS). Experi- 
ments were carried out in two greenhouses designed and realized 
in the Research and Technology Center of Energy (CRTEn) in 
Tunisia. We will present in Section 2, the potential of solar energy 
in Tunisia. In Section 3, we will describe the design of the 
experimental greenhouses, the experimental setup and measuring 
equipments. In Sections 4 and 5, a thermal analysis of the climate 
greenhouse and the energy balance of the SAHLSC are carried out. 
The thermal load leveling (TLL) are defined in Section 6. The 
experimental results are reported in the Section 7. Finally, we will 
present the economic analysis of the SAHLSC in Section 8. The 
main remarks of this work will be reported in the conclusion. 


2. Potential of solar energy in Tunisia 


The average amount of solar energy radiating in Africa is about 
6 kWh/m7/day, depending on the latitude of the country [34]. In 
the map of average global horizontal solar radiation (Fig. 1), we 
can see that there’s an important solar potential in North-Africa. 

Tunisia is one of the sunlight countries, it lies on the Medi- 
terranean coast of Africa; it is bordered by Algeria to the west and 
Libya to the southeast; with the following coordinates: Latitude 
36°43’N and Longitude 10°25’E. Tunisia's climate is temperate in 
the north, with mild rainy winters and hot summers. Tempera- 
tures in July and August can exceed 40 °C, while winters are mild 
with temperatures rarely exceeding above 20 °C [35]. 

The meteorological station on the site of Borj Cedria in Tunisia 
permits the measurement of the global sun flux on a horizontal plan, 
ambient temperature, wind speed and wind direction. Figs. 2 and 3 
give the variation of meteorological data of the average monthly global 
solar radiation and ambient temperatures in a period of 2009-2011. 
We can see that the global solar flux has the same trend of variation as 
the ambient temperature. It is observed that the peak of average 
monthly radiation is registered in July (225 kW/m?) with a highest 
monthly ambient temperature surrounding 30 °C. The least average 
monthly radiation is 65 kW/m? with a monthly ambient temperature 
around 13 °C in December. 
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Nomenclature 

A surface area (m?) 

Cp specific heat of air at constant pressure (J/kgK) 
da density of air (kg/m?) 

H heat loss coefficient of SAHLSC (W/m?°C) 
h heat transfer coefficient (W/m? °C) 

Ir total solar radiation (W/m?) 

L latent heat (J/kg) 

LAI leaf area index 

Loh characteristic length of greenhouse (m) 
L collector length (m) 

L2 collector width (m) 

Lz collector depth (m) 

l characteristic length of the leaf canopy (m) 
lea average length of the canopy (m) 

m mass (kg) 

m- mass flow rate (kg/s) 

Nu Nusselt number 

P(T) saturated water vapor pressure at temperature T (kPa) 
Pr Prandtl number 

Q heat rate (W) 

Ra Rayleigh number 

RE rate of air infiltration (m?/s) 

Re Reynolds number 

Ta aerodynamic resistance (s/m) 

Ts stomatal resistance (s/m) 

T temperature (K) 

t time (s) 

U velocity (m/s) 

V volume (m°) 

Z depth of soil (m) 


Greek symbols 


a absorbtivity for solar radiations 

at absorbtivity for thermal radiations 

ôb bottom insulation thickness (m) 

Ôe edge insulation thickness (m) 

A difference in time 

e emissivity 

y psychometric constant, 0.0667 kPa/K 

A thermal conductivity (W/m K) 

Ài thermal conductivity of insulation (W/m K) 
H PCM dynamic viscosity (Ns/m?) 

Hs PCM dynamic viscosity at T = Tpcy (Ns/m?) 
p reflectivity 


3. Materials and methods 
3.1. Experimental greenhouses description 


In this study, experiments were carried out in two small 
chapel-shaped greenhouses with same dimensions. The experi- 
mental greenhouses occupy a floor area equal to 14.8 m?, 3.7 m 
wide, 4m long and 3 m high at the center. The structure is all 
galvanized steel fixed to the ground with stones and concrete. 
The greenhouses are equipped with a centrifugal fan controlled by 
a differential thermostat. The fan operated when the temperature 
inside the greenhouse exceeds the optimal growth temperature of 
plant. A tomato crop (local variety) was planted in both green- 
houses and planted at 7 plants in the row and arranged in 4 rows 
with 40 cm between the rows and 35 cm between the plants. 


o Stefan—Boltzmann constant, 5.670 x 1078 W/m? K* 
T transmissivity 

Subscripts 

A absorbed 

a ambient 

av average 

G cover of greenhouse 

ca canopy 

b bottom 

ch charging 

dis discharging 

e edge 

fin_ch final of the charging process 


fin_dis final of the discharging process 
g glasses of SAHLSC 

H heat gain 

i inside greenhouse 
in inlet of SAHLSC 

inf infiltration 

L liquid phase 

los leakage, losses 

m melting 

max maximum 

min minimum 

o outsider greenhouse 


out outlet of SAHLSC 


PCM phase change material 
p packed-bed absorber 
S solid phase 

S soil 

sky sky 

st storage 

t top 

u useful 

v vertical section of the greenhouse 
Exposants 

A absorbed heat 

C convective heat 

Cd conductive heat 

L latent heat 

R radiation heat 


Physical characteristics of the greenhouse glass cover, canopy and 
soil are listed in Table 1. The reflectivity and transmissivity spectra of 
the canopy and cover layer recorded in the wavelength range (380- 
1400 nm) were plotted in Figs. 4 and 5. respectively. The reflectivity 
and transmissivity spectra were realized using LAMBDA950 UV-vis- 
NIR spectrometer equipped with an integrating sphere. 

The first greenhouse is named a Conventional Greenhouse (CG) 
where its walls and roofs are covered with a 3 mm of thick glass 
(Fig. 6). The second used a heating system and was named as an 
Insulated Greenhouse with Latent Heat System (IGLHS) (Fig. 7(a)). 
The south oriented wall and roof of the IGLHS are covered with a 
3 mm of thick glass. Sidewalls and the northern roof are built with 
0.4 m and 0.6 m of thick sandwich panels respectively. 

The heating system of the IGLHS is a new Solar Air Heater with 
Latent Heat Storage Collector (SAHLSC) using a packed bed of 
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spherical capsules as a latent heat storage system [36]. A new type 
of the SAHLSC was designed as a means to heat the interior 
environment of the greenhouses during the nighttime. Fig. 7 
(b) shows the SAHLSC manufactured inside the IGLHS. The length, 
the width and the total volume of the collector are 2 m, 1 m and 
0.28 m°, respectively. A 0.004m of thick transparent glass 
cover was placed 0.015 m apart from the packed bed absorber. 
A 0.05 m thick polyurethane insulation, with heat conductivity of 
0.028 Wm~'K~', is placed in the bottom and the edge of the 
collector. The packed bed absorber is the most important compo- 
nent of the solar air heater collector, formed of spherical capsules 
with a black coating and fixed with steel matrix. Capsules, 
confined with a PCM, have an outer diameter of 0.077 m and are 
blow molded from a blend of polyolefin with an average thickness 
of 0.002 m [36]. 


3.2. Experimental setup and measuring equipments 


Experimental greenhouses (east-west orientation), are con- 
structed and located in the Research and Technology Center of 
Energy in Borj Cédria, in Tunisia. 

During the sunshine time (charging process) a fraction of the total 
solar radiation received inside the greenhouse is absorbed by the 
black packed bed absorber surface of the SAHLSC. The absorbed 
thermal energy is stored as a sensible and a latent heat forms into 


o> a ie 


the collector. At night, a fan blows air across the PCM capsules 
and extracted the stored heat into the IGLHS, the discharging 
process done. 

K-type thermocouples inserted inside and outside the greenhouses 
for measuring the horizontal distributions of temperature, the vertical 
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Fig. 2. Monthly variation of the global solar radiation in Tunisia (2009-2011). 
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Fig. 1. Global solar radiation map of Africa [34]. 
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Fig. 3. Monthly variation of the ambient temperature in Tunisia (2009-2011). 


Table 1 
Cover, canopy and soil characteristic. 


Cover Canopy Soil 
Reflectivity for solar radiations 0.10 0.1 0:2 
Reflectivity for thermal radiations 0.10 0.55 0.16 
Emissivity 0.9 0.95 0.9 
Transmissivity for solar radiations 0.85 0 0 
Transmissivity for thermal radiations 0.88 0.4 0 


Reflectivity (%) 
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Fig. 4. Reflectivity spectra of canopy and cover. 


distributions and the temperatures at the top and the bottom of glass 
and sandwich panel. The temperature and relative humidity inside 
greenhouses (HMP155A, Campbell Scientific Inc) were measured in 
the center at a height of 1.5 m from the ground. The HMP155A sensor 
was instaled inside a 14-plate radiation shield houses and protects. Its 
louvered construction allows air to pass freely through the shield, 
thereby keeping the probe near ambient temperature [37]. Tempera- 
tures under the soil inside the greenhouse at 0, 0.25 and 0.5 m depth 
were measured using three PT-107 sensors. The average temperature 
of the canopy was measured using a non-contact means of measuring 
the surface temperature by an infrared temperature sensor IR120. The 
IR-SS Solar Shield is used to protect the IR120 sensor from direct solar 
radiation [37]. The global solar irradiation in the horizontal plane was 
measured with a Kipp and Zonen pyranometer and situated at 1.5m 
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Fig. 5. Transmissivity spectra of canopy and cover. 


Fig. 6. External view of conventional greenhouse. 


above the ground in the center of the IGLHS. Ambient temperature 
and relative humidity are measured by a HMP155A sensor situated 
at a height of 1.5 m above the greenhouse. Temperatures along the 
SAHLSC were measured using K-type thermocouples. Four sensors 
were inserted inside the capsules along with the absorber packed bed. 
Also four thermocouples are used to measure the absorber surface of 
one capsule, the glass, the inlet and the outlet air temperatures of the 
SAHLSC. All climatic parameters and measuring instruments were 
sample recorded every 10 min using a CR5000 data logger (Campbell 
Scientific Inc). 

The external global solar radiation, wind speed and direction 
are provided by the NRG weather station in the Research and 
Technology Center of Energy in Tunisia. This station is equipped 
with an acquisition system that record, every 10 min the average, 
the maximum, the minimum and the standard deviation values for 
each sensor. The weather station provided a wind speed measure- 
ment at 20 m. The wind velocity at the altitude of 2 m above the 
greenhouse was corrected by using the power law equation [38] 


(U1 /U2) = (hy /h2)? ad) 


where U and U, are the wind speeds at heights h; and hz 
respectively. # is the power law exponent and the subscripts (1 
and 2) correspond to different heights. 

Uncertainty analysis is needed to prove the accuracy of the 
experiments. The independent parameters measured in the 
experiments reported here to temperature, air velocity, and solar 
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Fig. 7. (a) External view of IGLHS and (b) SAHLSC inside greenhouse. 


radiation. To carry out these experiments, the sensitiveness of data 
acquisition system is about +0.001 °C, the measurement error 
is + 0.002 °C, the sensitiveness of the thermocouple is + 0.01 °C, 
the sensitiveness of the IR120 is + 0.0004 °C and sensitiveness of 
the PT107 is + 0.01 °C. The HMP155A errors are + 0.02 °C about 
temperature and +2% about humidity. An anemometer with 
+0.01 ms~! accuracy, and Kipp and Zonen pyranometer with 
+ 3% measurements uncertainties were used. The sensitiveness 
was obtained from a catalog of the instruments. 


4. Thermal analysis of greenhouse 


To describe the thermal behavior of the greenhouse, there are 
four components that play important roles in the thermal balance: 
the cover, the canopy, the soil and the inside air are analyzed. The 
following assumptions are made: 


— Radiation heat exchange between the greenhouse walls and 
roofs is neglected. 

— Storage capacity of the sandwich panel walls and roof material 
is neglected. 

— Absorbtivity and heat capacity of the enclosed air is neglected. 

— Conduction heat exchange between air and the sandwich panel 
walls and roof of the greenhouse is neglected. 


4.1. Energy balance of greenhouse cover 


The energy balance of the cover of greenhouse described by the 
following equation: 


Os ae OF ¢ Or. QS i+ QF y+ Qe + Qh. =0 (2) 


where Qa is the solar radiations absorbed by the glasses of the 
greenhouse; OF. ¿ is the solar radiations reflected by the canopy 
and absorbed by cover; Q4_, is the reflected solar radiations by 
the soil surface and absorbed by cover; QS_, is the convective 
heat transfer from external cover of the greenhouse to ambient 
air; QS_; is the convective heat transfer from internal cover of the 
greenhouse to inside air; OF spy is the thermal radiation 
exchanged between cover and sky; Q*_. is the thermal radiation 
emitted by soil surface and absorbed by cover and Q*®_. is the 
thermal radiation emitted by canopy and absorbed by cover. 

The energy balance of the cover Eq. (2) can be written as 


IrAcac+ ITTcacpeaAca i IptcacpsAs = ne oAc( Tp —To)— ie. Ac(Tc _ Tj) 


+ TeAc(T sky = Ta +AsCarcesTs4 +AcaFatc&cal oy =0 (3) 


Different relations have been used to calculate the external 
heat loss coefficient at the outside surface of greenhouse covers 
[39-42]. It depends on the greenhouse geometry under any 
environmental conditions. Mac Adams [40] suggested a relation, 


which is used in this study 
he, =5.7+3.8U, (4) 


where U, is the wind speed outside the greenhouse. 

The evaluation of internal convective heat transfer coefficient 
follows Seginer and Livne [42] which considers the combined 
effect of free and forced convection. 


se U,\ 1/2 
he; =1.52|T.—To|/2+5.2 (7) (5) 

Ua = RE/A, is the average speed of the air inside the green- 
house; RE and A, are the rate of air infiltration and surface area of 
the vertical section of the greenhouse respectively. 


4.2. Energy balance of canopy surface 


The energy balance of canopy surface inside the greenhouse is 
described by the following equation: 


Qla + OF ca se Qazi On =0 (6) 


where Q4, is the shortwave radiation absorbed by canopy; Q* _,,, is 
the thermal radiations emitted by cover and absorbed by canopy; Qe 
is the thermal radiation emitted by canopy; Q°,_; is the convective 
heat transfer from canopy to inside air and Q/,, is the latent heat 
transfer between the canopy and the inside air of greenhouse. 

The energy balance for canopy in the greenhouse Eq. (6) can be 


written as 


tclracaAca + ecoTsAcattca = cad T4 Aca = hS, E iAca(Tca —T;) 
daCpLAI(P*(Tea)— PCT) 4 
yY Ta+ Ts 


ca = 0 (7) 


The evaluation of the convective heat transfer coefficient between 
the inside air and the canopy, Stanghellini [43] proposed 
Cc pen daCp 
i 117a? 


(lealT ca — Ti] +207U2)"4 (8) 


where Usa is the air velocity through the canopy inside the 
greenhouse. 


(9) 


The saturated water vapor pressure P*(T;q) at the temperature 
of the canopy is calculated by Tetens equation (Eq. (10)), the 
following empirical formula relating accurately saturated water 
vapor pressure to temperature for the temperature range between 
0 and 60 °C [44,45] 


(10) 


P*(Tca) = 0.6108exp( 17.277 ca ) 


Tca + 237.3 
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The aerodynamic resistance, rg (sm~!) of the canopy, mainly 
depends on the aerodynamic regime that prevails in the green- 
house and the leaf length. It is expressed as [45] 

p2 
Ta = 220—, 11 
a ups (1) 
where | is the characteristic length of the leaf and U; is the mean 
interior air speeds. 

The stomatal resistance rs (sm~!) of the canopy is derived from 
a simple empirical relationship, Eq. (12), with global radiation 
(thermal and humidity dependences of greenhouse plants tran- 
spiration were neglected) [45,46] 


- ) (12) 


rs=200(14+ ages; 


4.3. Energy balance of soil surface 


The energy balance of soil surface inside the greenhouse is 
described by the following equation: 


QS -Qs -Qs -Qs +Q- =0 (13) 
where Q4 is the absorbed solar radiation by soil; QS_; is the 
convective heat transfer from soil surface to air inside the green- 
house; Q© is the conductive heat transfer into soil surface; QË is 
the thermal heat radiation emitted by soil and Q®_, is the thermal 
heat radiation emitted by greenhouse cover and absorbed by soil. 


The energy balance of soil in the greenhouse (Eq. (13)) can be 
written as 


A 
IptcasAs —hS_;As(Ts —T;) SATs Tz) — esa AAs +ecoT4Acas = 0 
(14) 


The convective heat transfer coefficients between the air inside 
the greenhouse and soil surface given by Seginer and Livres [42] 


U..\ 1/2 
Ki =152T -T1452 (7) (15) 
Lgh 


4.4. Energy balance of air inside greenhouse 


The energy balance of inside air of greenhouse, only the 
convective heat transfer is considered. It is assumed that no 
thermal or solar radiation is absorbed by the air. The energy 
balance is described by the following equation: 


QE tQ itO int =0 (16) 


where QS_; is the convective heat transfer from inside air to cover 
of the greenhouse, QS,_; is the convective heat transfer from 
inside air to canopy surface, Q°_; is the convective heat transfer 
from inside air to soil surface, Q&; is the heat loss due to 
ventilation or infiltration and Qq is the heat input into greenhouse. 

The sensible heat balance of the inside air of greenhouse 
(Eq. (16)) can be written as 


hE Ac(Te—Ti) + NG _ jAca(Tca— Ti) +h _;As(Ts-T) 
—dgCpRE(T;—To)+Qy = 0 (17) 


where Q,, is the heat gain from the solar air heater with latent 
storage collector into the IGLHS. 


5. Energy balance of the SAHLSC 


The energy analysis balance during both the charging and 
discharging phases of the SAHLSC is 


Qa = Qu +Qst + Qios (18) 


where Qa, Qu, Qst and Qis are the absorbed, useful, stored and 
lost energies, respectively. 

Based on Duffie and Beckman [47], the useful heat gain from 
the collector is 


Qu = MCp(Tout — Tin) (19) 


The radiation absorbed flux of the absorber packed-bed is 
defined as 


Qa = AplagTg)lr (20) 
The stored heat flux during the charging and discharging 
phases is given by 
Qen = [MpcmCp,s(Tm — Tini_ch.pcm) + MpcmL 
+mMpcmCp I(T fin_cn pcm — Tm)]/Aten (21) 


Qais = [MpcmCp.L(Tm — Tfin_dis pem) + MpcmL 
+mpcmCp,s(Tin_dis pcm — Tm)]/Atais 22) 
The lost heat flux is given by Eq. (23), His is the collector 
overall heat loss coefficient. The thermal energy is lost from the 
collector to the surroundings by conduction, convection and 
infrared radiation. H),, is equal to the sum of energy loss through 


the top (H+), bottom (Hp) and edges (He) of the collectors given 
below [47] 


Qios = HosAc(Tp — Tj) (23) 


Hios = Ht + Hp + He (24) 
The top loss coefficient from the collector to the ambient is 
Hr=(1/(hp_g thp_,)+1/he_,)7' (25) 


The flow is laminar during the charging phase; and the 
appropriate correlation is given by Churchill [48] 


Nuch = 2+ (0.589 Ra'/*)/(1 +(0.469)/Pr)?/1°)*/° 
Pr > 0.7 
Ra <10!! (26) 


Then the flow is turbulent during the discharging phase; the 
appropriate correlation is given by Whitaker [44] 


Nugis = 2+(0.4 Re!/2 +0.06 Re?) Pr°4(u/p,)'/4 
0.71 < Pr < 380 
3.5< Re < 7.6 x 104 (27) 


1 <(u/ps) <3.2 


hp_ g =Nu(a/t) (28) 
hpg = (a(Tp+T2(Tp+Tz))/(1/ep) + (1 /eg)— 1) (29) 
hga = £g0 (T2 + Tay (Tz +T sty) (30) 


In addition the energy loss through the bottom (H;) and the 
edges (He) is as follows: 


Hp = 1i /Sp (31) 


He = (Ly + 1p) 134; /(Li L260) (32) 


6. Thermal load leveling (TLL) 


The thermal load leveling (TLL) is a relative index to represent 
the fluctuation of temperature inside the greenhouse. The less the 
fluctuations, the better is the environment for plants inside the 
greenhouse. For thermal heating of the greenhouse with minimum 
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fluctuations the TLL should has a minimum value [49]. 


T; ier 
TLL = imax i, min (33 
T;, max + T;, min ) 


7. Results and discussions 
7.1. Comparative study of greenhouses climate 


In the CG and IGLHS, tomato plants are transplanted on the 
23th of January and the harvesting started on the 15th of April 
2013. The major environmental factors that affect the growth and 
the precocity of the production of greenhouse plants are carbon 
dioxide, light, humidity and temperature. The tomato plants have 
two optimum temperatures, one during the day; which varied 
between 20 and 30 °C, and the other nocturnal temperature which 
is the most crucial temperature, varied between 12 and 18 °C. 

During the growing period, the February has the least average 
monthly radiation is 85 kW/m? with a monthly ambient temperature 
around 12.5 °C (Figs. 2 and 3). An experimental study has been carried 
out during the month of February 2013 which is the winter season. 

Fig. 8 presented the instantaneous registered air temperature 
inside CG and IGLHS (1st-28th February, 2013). The fluctuation of 
air temperature inside the two greenhouses is closely attributed to 
weather conditions such as solar radiation and outside tempera- 
ture. Under the experimental conditions, the nocturnal variation of 
temperatures inside greenhouses varied between 2 and 5 °C. The 
method of utilizing the new solar air heater with latent storage 
collector as a heating system inside the greenhouse shows mean- 
ingful value in February. 

An experimental study has been carried out on typical selected 
days of the system operation from 23th to 27th February, 2013 in 
winter season. 

The external global solar radiation and wind speed for the 5- 
day period in the region of the greenhouses are shown in Fig. 9. 
The maximum of the global solar irradiation intensity in the 
horizontal plane varied between 300 and 650 Wm~?. The two 
first days of this period were very windy, air speed varied between 
3 and 8 ms~'. The average speed of other days is of 2 ms~!. 

The daily variation of the ambient temperature, air tempera- 
ture inside CG and air temperature inside IGLHS are presented in 
Fig. 10 (for the same 5-day period). The air temperature inside the 
greenhouses followed the same trend as the ambient temperature 
and exceeds the daily tomato optimum temperature under severe 
winter weather. Consequently, it is necessary to improve the air 
conditioning of this greenhouse. A new solar air heater collector 
with a latent heat storage system is integrated inside the IGLHS. 
This collector is used to store the excess of the transmitted solar 
irradiation in the packed bed absorber and to provide it at night. 
The discharging process of the SAHLSC starts at 21 h (local time). 
A fan was used to blow the air at a fixed speed equal to 1m/s. 

The 23th and 27th of February, 2013 is characterized by clear 
sky with maximum ambient temperature equal to 21°C. We 
observed that the diurnal temperature attained 39°C and 41 °C 
respectively into CG and into IGLHS. The nocturnal variation of 
temperature inside the IGLHS exceeds the temperature inside CG 
with a difference of 5 °C from 21 h (local time) until the morning 
after (Fig. 10). On 27th of February, 2013 the SAHLSC maintained 
the air temperature inside IGLHS constant, around 12 °C, along the 
night, due to the fraction of solar radiations stored in the packed 
bed absorber. However, the temperature inside CG decreases 
progressively to the ambient temperature after 2 h from sunset. 
The 24th February, 2013 was characterized by severe solar radia- 
tion fluctuations, the ambient air temperature does not exceed 
11 °C and wind velocity varied from 4 to 8 m/s. We noted that the 
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Fig. 8. Air temperature inside CG and air temperature inside IGLHS as a function of 
days. (1st-28th February, 2013). 
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Fig. 9. Global solar radiation and wind velocity as a function of days (23th-27th 


February, 2013). 
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Fig. 10. Ambient temperature, air temperature inside CG and air temperature 
inside IGLHS as a function of days (23th-27th February, 2013). 


temperature into the IGLHS exceeds the temperature inside CG 
with a difference of 3 °C all the night. 

Fig. 11 presents the effect of the SAHLSC on relative humidity 
inside the greenhouse. At night, the IGLHS permitted to maintain 
the air relative humidity to 10-20% lower than the CG. The SAHLSC 
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Fig. 11. Humidity inside CG and IGLHS as a function of days (26th-27th 
February, 2013). 
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Fig. 12. Useful, absorbed and stored heat rates as a function of day times inside the 
IGLHS (27th February, 2013). 


creates a passive dehumidification process at night time due to the 
increase of the air temperature inside the IGLHS. 


7.2. Performance of the SAHLSC inside IGLHS 


Variations of useful, absorbed and stored energy are presented 
in Fig. 12 during 27th of February, 2013. Therefore, to evaluate the 
stored energy, the energy balance equation (Eq. (18)) has been 
used. A great part of the absorbed solar heat was stored inside the 
PCM. It is observed that during the initial period of charging, the 
instantaneous heat stored increases with insolation and towards a 
maximum value of 0.7 kW at 13:00, we noted that at this time the 
absorbed heat is 1.3 kW, 53% of the solar energy is stored into 
the collector. As the discharging process proceeds, the PCM 
starts solidifying and the used heat is uniform for a longer period. 
The uniform value of the used heat is between 550 and 300 W all 
the night. 

Energy balance equation (Eqs. (2), (6), (13), (17) and (18)) have 
been used for calculating requirements, excess and recovered 
thermal energy (MJ) inside IGLHS. Fig. 13 presented the daily 
heating requirements and the nighttime heating recovered from 
SAHLSC, during the 5-days period (23th-27th February, 2013). The 
stored energy in the solar air heater with thermal storage is 
13.7 MJ on 23th February, 2013. The contribution of stored energy 
in the SAHLSC to the nighttime heating requirements was 31%, 27% 
and 26% respectively on 23th, 26th and 27th February, 2013 under 
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Fig. 13. Nighttime heating recovered from SAHLSC and IGLHS heating require- 
ments as a function of days (26th-27th February, 2013). 
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Fig. 14. Thermal load leveling of CG and IGLHS as a function of days (26th-27th 
February, 2013). 


best ambient conditions. On 24th and 25th February, 2013, the 
recovered and stored heat is 6% and 11% with requirements heat 
exceeds 65 MJ. 

The performance of IGLHS with a SAHLSC heating system has 
been evaluated in terms of thermal load leveling, TLL, using Eq. (3). 
The TLL is used to quantify the fluctuations of temperature inside 
the greenhouse. In winter, TLL should have lower value by incor- 
porating a heating system due to the increase of (Tj, max +Ti, min) aS 
well as decrease of (Ti, max—Ti, min) as Compared to TLL without 
heating system. The results of the daily variation of thermal load 
leveling for CG and IGLHS greenhouse have been shown in Fig. 14. 
It can be seen that the value of TLL is maximum for CG greenhouse 
and it is reduced to 25% for IGLHS on 27th February, 2013. The 
lower values of thermal load leveling indicate the decrease in the 
fluctuations of greenhouse air and thereby, there occurs an 
improvement of desired environment for plants in the greenhouse. 


8. Economic analysis 


The economic analysis (SAHLSC cost, energy cost and payback 
period in comparison with the expected service life cycle of the 
installation) is done without taking into account the interest and 
the inflation, we suppose that all the cost is paid without credit 
payment [50]. The life cycle of the SAHLSC is considered 20 years 
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and its thermal performance degradation is assumed to be 1% 
per year. 

In order to evaluate the economic analysis, the SAHLSC system 
is compared to non-solar alternative heating system (fuel boiler). 
In this study, two systems are considered: 


— Hybrid system: the SAHLSC and the fuel boiler. The fuel boiler 
is the auxiliary system meeting the energy needs. 

— Conventional system: the most heating system used for green- 
houses conditioning is the fuel boiler. 


8.1. The total cost system 


The total cost of the greenhouse heating system includes the 
investment cost and the energy consumption cost. The investment 
cost includes the system cost and the operation and maintenance 
cost [51,52]. 

The total cost of the conventional system is defined by adding 
the initial investment cost of the fuel boiler and the fuel con- 
sumption. For the hybrid system the total cost is the sum of all the 
initial investment cost (SAHLSC+fuel boiler) and the fuel con- 
sumption without the supplied heat cost by the SAHLSC. 


— The system cost is defined by the equipment, the realization 
and the installation costs. The hybrid system cost includes the 
fuel boiler and SAHLSC costs. Since the SAHLSC is not yet 
commercialized, there is no official price for this system. Thus, 
the capital cost of the solar collector was roughly determined 
as the costs spent to develop the experimental system. 

— The operation and maintenance cost refers to heat production 
coupled with full load hours. The electricity consumption for 
operating the fan was considered as operating costs. The 
maintenance cost of the SAHLSC is considered negligible. 


8.2. Energy consumption for the greenhouse heating 


The sizing of the greenhouse heating system is based on the 
nocturnal setpoint temperature (12 °C). This temperature is sui- 
table for tomato cultivation implemented in the experimental 
greenhouse. The necessities of the energy to maintain the noctur- 
nal air temperature under the greenhouse at the setpoint value is 
equal to the energy losses. The total energy losses are determined 
by adding the energy losses of the greenhouse components 
(cover, ground, infiltration, and canopy). The Tunisian greenhouse 
heating period is from November to March, equivalent to 1500h 
of heating. In this period, the fuel boiler in the hybrid system, 
functions only 600h. Therefore it ensures an economy of 
1800 kWh/year. 


8.3. The payback duration 


The aim of the economical analysis is to find out the payback 
period of the SAHLSC. The payback duration indicates the time 
that the SAHLSC generates the amount of energy equivalent to the 
energy provided during the life cycle [53]. 

The cost of each heating systems (hybrid and conventional 
system) and the payback time of the SHALSC is mentioned in 
Table 2. As a result, the payback period of solar air heater with 
latent storage is found as 5 years. 


8.4. Environmental impact of the SAHLSC 


From 2005 to 2008, the greenhouse areas of Tunisia increased 
from 4600 ha to 9000 ha [54], which consume 3 x 10° t of fuel 


Table 2 
Comparison of the cost and the payback times between SAHLSC system and the 
fuel boiler. 


Hybrid system Conventional system 


(SAHLSC + fuel boiler) (fuel boiler) 

System cost ($) 3600 1800 
Operation cost ($) 4.75 116.55 
Maintenance cost ($) 150 150 
Annual energy 1200 3000 

consumption (kWh) 
Gain of energy (kWh) 1800 - 
Energy cost ($) 259.2 651.25 
Total cost 4014 2717.8 
Payback of SAHLSC 5 - 

(year) 


during 5 months of heating period. The maximum of the CO, 
content is about 16% of the fuel consumption [55], therefore the 
conventional heating system causes the emission of 3.12 x 10° t of 
CO, which is equal to 2.8 x 107? t per person [56]. The reduction 
of carbon dioxide emissions using the SAHLSC may be a promising 
solution. Indeed, the SAHLSC ensure the greenhouse conditioning 
during 3 months, it can reduce the CO, emission at 1.68 x 10-7 t 
per person which presents 1% of the total of the CO, emissions. 


9. Conclusion 


An experimental study was conducted to evaluate the night- 
time recovered heat of the SAHLSC in the east-west oriented 
greenhouse. During winter months, the maximum of solar radia- 
tion falls on the south wall. Therefore, a SAHLSC is proposed in this 
study for storing the excess of the transmitted solar irradiation in 
the packed bed absorber and to provide it at night. This stored 
energy is extracted to the greenhouse by convection and radiation 
heat. The performance of the solar system was evaluated in terms 
of thermal load leveling and nighttime recovered heating for 
a 5-day period of February, 2013 in Tunisia. As a result of this 
system, the amount of the nighttime recovered heat of this system 
attains 31% of the total requirements of heating. The relative 
humidity was found to be an average 10-20% lower at night time 
inside the heated greenhouse. The nocturnal variation of tempera- 
ture inside the IGLHS exceeds the temperature inside the CG by 
5 °C. The solar air heater remains a uniform useful heat during the 
nighttime. The uniform value of the used heat is between 550 and 
300 W all the night. This useful heat was not affected by the severe 
global solar radiation fluctuation during the charging mode. The 
evaluation of the economic analysis of the SAHLSC assessed that 
the payback period was approximately 5 years if the system 
operates only 3 months per year. The environmental impact from 
reducing carbon emissions is about 1% of the total Tunisian 
greenhouses emissions. 
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